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Abstract 

Several classes of systems of evolution equations with one or two vector 
unknowns are considered. We investigate also systems with one vector and 
one scalar unknown. For these classes all equations having the simplest higher 
symmetry are listed. 

1 Introduction. 

The symmetry approach to the classification of integrable evolution equations (see 
[H El 131 Hj) is proven to be the most efficient integrabihty test for 1+1-dimensional 
nonhnear PDEs. 

In this paper we extend the simplest version of this approach |5l |6] to the case 
of so called vector evolution equations (see [7]). As an example let us consider the 
following vector equation 

Ut = U,xx + {U,U)U,, (1) 



where U{t,x) is a A^-component vector and (■,■) stands for the standard scalar 
product. If = 1 this formula reduces to the well-known modified Korteweg- 
de Vries (mKdV) equation ut = Uxxx + u^Ux- Equation ([1]) belongs to a wide class 
[3 13 m [TOl [11], [12] of integrable multi-component evolution systems related to Jordan 
algebras and Jordan triple systems (see Appendix 1 for definitions). 

S.I. Svinolupov was the first who discovered a deep relationship between some 
kinds of non-associative algebraic structures and integrable multi-component evo- 
lution systems. In the next two pages we formulate his results that have relations 
with the subject of our paper. 

One of his results states that for any Jordan triple system {X, Y, Z} defined on 
a vector space J the following equation 

Ut = Uxx. + {U,U,Ux}, (2) 

where U{x, t) is an unknown function with values in J is integrable (and, in par- 
ticular, has higher symmetries). To derive some concrete vector integrable systems 
from this general result we remind that there exist two different "vector" Jordan 
triple systems: 

{X, Y, Z} = (X, Y)Z+ {Y, Z)X - (X, Z)Y (3) 

and 

{X,Y,Z} = {X,Y)Z+{Y,Z)X. (4) 

Here X, Y and Z are arbitrary X-dimensional vectors. It is clear that ([3]) gives rise 
(up to an inessential scaling parameter) to ([1]) while (jlj) yields a different integrable 
vector generalisation of the mKdV equation 

Ut = Uxx. + {U, U) Ux + (f/, Ux) U. (5) 

The general formula ([2]) produces not only one-component vector equations of 
mKdV-type. Let us consider the following Jordan triple system 




n{Y\z^) + {Y\ z^)) X' +{{x\ z') - {x\ z^)) +((xi, yi) + {x\ y^)) z^\ 

\{{Y\ Z') + {Y\ Z^)) X^ +{{X\ Z^) - {X\ Z^)) Y^ +((Xi, Y') + {X\ Y^)) Z^j 

(6) 

where X^, Y^ , Z^ are vectors of length Ni and X^, y^, Z^ are vectors of length N2. 
This Jordan triple system was not mentioned in the original papers by Svinolupov. 
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The corresponding system 

= Vl^^ + [{U\ U') + {U\ U')) Ul + 2{U\ Ul) 

U? = + {{U\ U') + {U\ U^)) Ul + 2{U\ Ul) U^ 

with respect to vectors U^ and U"^ of lengths A^i and N2 was presented in 
although in non-explicit form this system was probably found first in [13j. Notice 
that if A''! = A^2 = 1 then ([6]) is a fully decoupled system in variables f/^ = f/^ + U"^ 
and t/^ = U^ — U"^. For A^i = A^2 > 1 it is impossible to split the system by this 
change of variables. 

One more general result by Svinolupov is related to a mult i- component gener- 
alisation of the nonlinear Schrodinger equation (NLS) which can be written as a 
system of two equations 

Ut = u^x + u^v, Vt = -v^x - v^u. (7) 

It was shown in [9] that for any Jordan triple system the following system of equa- 
tions 

(Ut= Uxx + {U,V,U}, 
\Vt = -Vxx-{V,U,V} 

is integrable. In particular, it has higher symmetries of any order n G N. For the 
vector triple systems ([3]) and (jl]) the general formula ([8]) gives rise to two different 
integrable vector generalizations of ([7]): 

(Ut= Uxx + 2{U,V)U-{U,U)V, 
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I Vt = -Vxx - 2{U, V)V+ {V, V)U, 

and 

(Ut= Uxx + {U,V)U, 

{vt = -Vxx-{u,v)v. ^° 

Both systems are well known (see [IH[I5]). 

Finally, Jordan triple system ([6]) yields the following 4- vector system of NLS-type 



ul = Ul, + 2{{U\V') + {U\ V^)) f/i + {-{U\ U') + {U\ U')) V\ 
Ul = Ul, + 2{{U\V') + {U\ V^)) f/2 + {{U\ f/i) - {U\ U^)) V\ 
VI = -v,l - 2i{U\V') + {U\ V^)) V^ - {-{V\ V^) + {V\ V^)) u\ 
V^ = -VI - 2{(JJ\V^) + {U\ V^)) V^ - {{V\ V^) - {V\ V^)) u\ 



(11) 



where vectors and have length A'^i whereas U"^ and V"^ are vectors of length 
A^2- It is easy to see that the reduction U'^ = V"^ = reduces the system to ([2]). If 
Ni = N2 then the reduction f/^ = — f/^, V"^ = —V^ becomes possible and results in 
system (|TOl) . 

For all examples presented above, the right-hand side of the equation is a homo- 
geneous differential polynomial under a suitable weighting scheme. The differential 
equation 

Ut = f{u,U^, . . . ,Un-l,Un), ^ ^ ^^'^'^ 

is said to be X- homogeneous of weight /i if it admits the one-parameter group of 
scaling symmetries 

(x, t, u) — > (a~^x, a~^t, a/^u). 

For A^-component systems with unknowns -u^, -u^ the corresponding scaling group 
has a similar form 

In this paper we consider the case Ai = ■ ■ ■ = Aat only. 

It was proven in [16] that in the scalar case a A-homogeneous polynomial equation 
with A > may possess a homogeneous polynomial higher symmetry only if 

• Case 1: A = 2; 

• Case 2: A = 1; 

• Case 3: ^ = |- 

For example, the KdV equation ut = Uxxx + uu^ is homogeneous of weight 3 for 
A = 2, the mKdV equation ut = u^xx + u^Ux has the weight 3 for A = 1 and for the 
Ibragimov-Shabat equation 

Ut = Uxxx + "iu^Uxx + ^Uul. + "iu^Ux (13) 

the weight is 3 and A = |. 

According to jlTl [18] equations like f|T3|) do not exist in the matrix case. Thus 
the possible values of A for matrix one-component equations possessing higher sym- 
metries are 2 or 1. 

In the paper [21] two-component systems of the form 

Ut= kiUxx + P{ux,Vx,u,v), 

Vt= k2Vxx + QiUx,Vx,U,v), 
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where ki ^ k2, ki ^ 0, have been considered. The authors assumed that P and Q 
are polynomials such that the system is (Ai, A2)-homogeneous, where A2 > Ai > 
and A2 — Ai is not a natural number. Under these assumptions it was proven that 
if the system possesses higher symmetries, then 

• if both P and Q have no quadratic terms, then k2 = —ki, 

• if Ai = A2 = A, then A G {2, 1, i}. 

Also any possible pairs of Ai and A2 were listed. 

For integrable systems with more than two unknowns there are no rigorous state- 
ments describing all possible types of homogeneous polynomial systems but for all 
known examples with Ai = ■ ■ ■ = Aat = A > 0, A is equal to 2, 1 or ^ just as in the 
scalar and two-component cases. 

In this paper we systematically investigate one and two component homogeneous 
vector equations having higher symmetries under the assumption AG{2,l,i}. It is 
remarkable that the number of different integrable cases for vector equations turns 
out to be greater than in the matrix case. 

Methodologically our approach is not new. But the continuous increase of com- 
puter speed and development of efficient computer algebra algorithms, especially 
for solving overdetermined algebraic systems allows to tackle problems which were 
unsolvable so far. 

Nevertheless, even modern technology cannot deal with multicomponent sys- 
tems, if we use a straightforward component formalism. A corner stone of our 
approach in this paper is a special component-less representation of vectors and 
vector equations. 

Despite of such customised data-representations the paper is the result of long 
and intensive computations with the help of the computer algebra package Crack 
fi9[ [20] . Although originally designed to solve overdetermined PDE-systems, it got 
enhanced in recent years to solve algebraic and especially bi-linear systems efficiently 

We find a number of new vector models supposed to be integrable. The best way 
to prove the integrability of these models is to find Lax representations for them. 
We are planning to do that in a separate paper. 

2 Integrable equations with one vector unknown 

In this section we consider vector equations of the form 

Ut = fnUn + fn-1 + " " " + /l U, + foU, = ^ (14) 
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with a single unknown vector-valued function U of dimension A^. Here are scalar 
functions of scalar products {Ui, Uj), < i < j < n. We shall call equations of this 
form isotropic. It is clear that the isotropic equations are invariant with respect to 
any orthogonal transformation of the vector U i.e. the equation has the orthogonal 
group 0{N) as a group of point symmetries. 

We shall consider equations f[T^ that are integrable for any dimension A^. In 
addition, we assume that the coefficients fi do not depend on A^. In virtue of the 
arbitrariness of A^, all different scalar products {Ui, Uj) with i < j can be regarded 
as independent variables. Their functional independence is a crucial requirement in 
all our computations. If A^ would be fixed, we could not assume that. For instance, if 
A^ = 3, then the determinant of the matrix with entries aij = {U, Uj), i,j = 1, . . . , 4 
identically equals to zero. 

The signature of the scalar product is unessential for us. Furthermore, the as- 
sumptions that the vector space is finite-dimensional and the constant field is M are 
also unimportant. For instance, U could be a function of t,x and y and the scalar 
product be 

/oo 
U{t,x,y)V{t,x,y) dy. 

In this way, our formulas and statements are valid also for this particular sort of 
1 + 2-dimensional non-local equations. 

We suggest that the coefficients fi are polynomials in {U, Uj) such that (HM is 
homogeneous with A = 2 (Case 1), A = 1 (Case 2) or A = | (Case 3). Given the 
weight of f|T^ it is easy to find the most general form of such an equation (we have 
a special symbolic code for that). Everywhere we assume that /„ = const 7^ 0. 

For example, the general form for an equation of weight 3 with A = 1 is given by 

Ut = a,U,^^ + a^iU, U)U, + a^iU, U,)U, (15) 

where a, are arbitrary constants, ai 7^ 0. Using a scaling of t we can normalize ai 
to 1. 

To describe all integrable cases we assume (cf. with the scalar case [51 \3\ \T7[\W\ ) 
that (1151) has a higher homogeneous symmetry of weight 5, whose general form is 

Ur = Uxxxxx + hi{U, U) Uxxx + h2{U, Ux) Uxx + 

MU, + &4(f/x., U.,) + h,{U, U)^) U, + (16) 
ih{U, Uxxx) + br{Ux, Uxx) + bs{U, U){U, Ux)) U 

Compatibility conditions (t/t)r = (f^r)t of (IT^ and (ITBl) give rise to a system of 26 
bilinear algebraic equations with respect to unknowns 02, 03, 61, . . . , fcg- Solving this 
system with the help of the symbolic computer program Crack we finally obtain 
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that the system has a solution only if 03 = or 02 = — 203. These two possibilities 
lead to integrable equations ([T]) and 

We considered Cases 1,2 and 3 corresponding to different weighting schemes. For 
each case we have performed the similar calculations to find all equations of second 
order having a symmetry of third order, equations of third order having a symmetry 
of fifth order, and equations of fifth order possessing a symmetry of seventh order. 
The complexity of the algebraic system increases drastically with an increase of the 
total weight (differential order of the equation and symmetry) and with decreasing 
A. For details see Appendix 2. 

It turns out that apart from equations ([1]) and ([5]), there exists only one more 
integrable equation: 

Ut = U,,, + 3{U, U)U,, + 6{U, U,)U, + 3{U, U)'U, + 3{U,, U,)U (17) 

with A = |. This vector analogue of the Ibragimov-Shabat equation (fT3l) was already 
found in [22] . 

All equations and (ITTIl belong to infinite commutative hierarchies of 

evolution equations such that any equation from the hierarchy is a symmetry for 
all others. These three hierarchies are of the same (3,5)-type. This notation means 
that the simplest non-trivial equation in the hierarchy is of order three and the next 
is an equation of order 5. 

We see, that there are no hierarchies of type (2,3) and (5,7) in the vector case 
with one unknown vector in contrast with the scalar case, where such hierarchies 
exist (Burgers, Kaup-Kupershmidt, Kupershmidt and Sawada-Kotera equations). 

Combining our results with the approach by Sanders- Wang [161 IH] one can 
prove that any polynomial homogeneous equation of the form (IT^ with A > and 
fn = const 7^ which has at least one higher symmetry, belongs to one of the above 
three hierarchies [23] . 



3 Integrable equations with two vector unknowns 

In this section we classify integrable vector NLS-type systems. More precisely, we 
consider systems of the form 

Ut = U^x +PiUx+ P2K + P3U + PaV, 

(18) 

Vt = -Vxx + P^Ux + PgVx + PjU + psV, 

where U and V are vectors and the coefficients Pi are A-homogeneous polynomials 
depending on all possible scalar products of vectors [/, V, ^Tj, V^. By analogy with 
the scalar case A is supposed to be 2,1, or |. To derive the integrable cases we 
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assume that (fTSI) . just as in the scalar case (see [2]) possesses a symmetry of the 
form 

(Ur = U^xx + qiUxx + q2Vxx + q^JJx + qiVx + q^U + ge^, 

< (19) 

I K- = Vxxx + gyf^x + qsVxx + ggf/x + Qio^x + qiiU + gi2V, 

where the coefficients are A-homogeneous polynomials of all possible scalar prod- 
ucts of U^V^Ux^Vx^Uxx^yxx- We call such systems integrable. The classification 
result is the following. 

It is easy to see that nonlinear systems f|T8|) corresponding to the case A = 2 do 
not exist. 

For A = 1 there exist only two integrable systems (|TOl) and ([9]) (up to a scaling 
of f/, V, X and t). 

In the case A = | the complete list of integrable cases (up to the scaling and the 
involution U ^ V, t ^ —t) looks as follows: 

{Ut= Uxx + 2a{U,V)Ux + 2a{U,Vx)U -a(3{U,VyU, 
\Vt = -Vxx + 2(3{U, V)Vx + 2(3{V, Ux)V + a(3{U, VfV 



Ut = Uxx + 2a{U, V)Ux + 2^{U, Vx)U + /3(a - 2/3)(f/, VfU, 
Vt = -Vxx + 2a{U, V)Vx + 2(3{V, Ux)V - (3{a - 2(3){U, VfV 



(21) 



Ut = Uxx + 2a{U, V)Ux + 2(3{U, Vx)U + 2{(3 - a){V, Ux)U 
-a(]{U,V)^U, 

Vt = -Vxx + 2a{U, V)Vx + 2a(1/, Ux)V + VfV 



(22) 



f f/t = Uxx + 2a([/, + 2a(f/, + 2/3(y, Ux)U 
-a{a - P){U,VyU, 

Vt = -Vxx + 2a{U, V)Vx + 2a{V, Ux)V + 2(3{U, Vx)V 
+a{a- (3){U,V)^V 



(23) 



f Ut = Uxx + 4a(f/, V)Ux + 2{a - /?)([/, [/)V; + 4/?([/, Vx)U 
+Af3{a - 2/3) (?7, V')^^/ - 2(3{a - (3){U, U){V, V)U 
-AP{a-P){U,U){U,V)V, 

Vt = -Vxx + MU, V)Vx + 2{a - P){V, V)Ux + A(3{V, Ux)V 
-A(3{a - 2/3) (f/, V)'^V + 2[3{a - f3){U, U){V, V)V 
+APia-P){V,V){U,V)U 



(24) 
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' Ut = U^^ + Aa{U, V)U^ - 2I3{U, U)V^ + Aa{U, V^)U 

+4/3(y, U.,)U - 4/3([/, U,)V + 6l3{a - P){U, U){V, V)U 
-4a{a - p){U, VfU - 4/5(a - 0){U, U){U, V)V, 

' Vt = -Kx + MU, V)V, - 2P{V, V)U, + 4a{V, U,)V ^^^^ 
+4f3{U, V,)V - 4(3{V, V,)U - 6f3{a - P){U, U){V, V)V 
+4a{a - P){U, V)^V + 4(3{a - p){V, V){U, V)U 

All equations in the list contain arbitrary constants a and {3. It is easy to see 
that if a is not equal to zero, then it can be reduced to 1 via scalings of t, x, U and 
V . Thus the essential parameter is the ratio of (3 and a. We choose the above form 
with the two arbitrary constants to avoid considering the case a = separately. 

Remark 1. In the paper equations fl^T]) and fl2^ were found and Lax 
representations for them were presented. The other equations seem to be new. 

Remark 2. In the paper [27] a list of matrix integrable equations have been 
obtained. This list contains all integrable polynomial homogeneous matrix equations 
with A = Some particular cases of our vector systems (l20l) - (!25l) can be derived 
from the matrix list as reductions. For example, for all matrix systems we can 
perform the following reduction 
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to get the corresponding systems with four vector unknowns V, X, Y. Sometimes 
we can decrease the number of vector unknowns to two by further reductions and 
obtain some systems from the list (!20!) - (l25l) . But since non of the matrix equations 
from [27] have arbitrary parameters and our equations do, we probably can con- 
struct only particular equations from the list by reductions of the matrix 
equations. 



4 Integrable equations with one scalar and one 
vector unknown 

Some vector analogues of the KdV equation can be constructed using a different 
general result of Svinolupov [8J: for any Jordan algebra (see the definition in the 
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Appendix 1) with multiplication o the following "Jordan" KdV equation 

Ut = f/xxx + UoU^ (26) 

has higher symmetries. Since for any Jordan triple system J = {X, Y, Z} and given 
fixed vector C the following multiplication 

XoY = {X,C,Y} 

defines a Jordan algebra, the formulas ([3]) and (j4]) produce two vector KdV equations 

Ut = t/... + {U, C) + (C, U^) U - {U, U^) C (27) 

and 

Ut = U^.. + {U, C) U, + (C, U,) U (28) 

The first had been obtained in |7j and the second was considered in 
The Jordan triple system gives rise to 

f Ul = U'^^^+ {{U\C') + {U\C'))U'^ + 
{{C\Ul) + {C\Ul))U'+ 
{-{U\Ul) + {U\Ul))C\ 

u! = ul^+ {{u\c') + {u^c^))u'^+ 

{{C\Ul) + {C\Ul))U^+ 
{{U\Ul)-{U\Ul))C\ 

Note that it is impossible to construct any vector KdV equation (with a term 
quadratic in U) without using constant vectors. Indeed, the simplest nonlinear 
vector term which can be generated with the help of the scalar product from only 
[/, Ux, Uxx is of third degree. 

For non- isotropic equations like (1251) one can use the orthogonal group 0{N) to 
simplify the vector C and bring it, for example, to the form 

C= (1,0,0,...,0). (30) 

With vector C fixed, the equation admits the symmetry group 0{N — 1). 
The component form of fl28|) with C given by fl30|) is the following 

Ut = ulxx + U^ul + M^M^ 
N Af I 1 TV I 1 TV 



(29) 
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where U = {u^, . . . , u'^). We see that the first component satisfies the standard 
scalar KdV equation and the remaining equations are hnear if the function is 
aheady found. We call such vector systems triangular. 

In contrast, fl271) cannot be splited by any orthogonal transformation, i.e. it is 
non-triangular. Note that if C is of the form fl30l) then (|27|) can be rewritten as 

Ut = Uxxx + UUx - {U,Ux), 

Ut = Uxxx + uUx + UxU, 

where u = and U = (m^, . . . , -u^). 

System (!3T|) gives us an example of an integrable systems with one vector and 
one scalar unknown. Surprisingly, for such systems the list of integrable models 
seems to be much richer than in the case of isotropic equations. In the following a 
lower case u stands for the scalar variable and the vector variable is denoted by a 
capital U. 

4.1 The case A = 2 

In this paper we consider the simplest case A = 2. The generic form for a second 
order homogeneous system of this type reads as follows 

Ut = ai Uxx + 02 + as (J7, [/) , 

(32) 

Ut = 04 Uxx + a^uU, 

where Oj are constants. We assume that at least one of the constants ai or 04 is 
non-zero. It means that the system is non-degenerate. As in the scalar case, the 
classification result is negative. It turns out that system fl52]) has a symmetry of 
third order only if 02 = as = 05 = 0. 

For third order systems with A = 2 the result is much more interesting. The 
generic form of such systems is given by 

Ut = ai Uxxx + 0-2 UUx + a3{U, Ux), 

(33) 

Ut = 04 Uxxx + a5uUx + ae Ux U. 

These systems can be regarded as natural vector generalizations of the KdV equa- 
tion. We assume that 1) the system is of third order, i.e. at least one of the 
constants ai or 04 is non-zero; 2) our system is non-triangular, i.e. as 7^ and one 
of 05 or Og is non-zero. 

Theorem. Suppose a non-triangular system (1331) possesses a symmetry of fifth 
order. Then (up to a scaling of t,x,u and U) this system coincides with (1311) . or 
belongs to the following list 
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Uxxx + 3uux + 3{U, Ux) 



(34) 




Uxxx + UU^ + ([/, f4), 



{U, u,) 



2a 



XXX 



(35) 



(36) 



System fl34|) has been considered by Kupershmidt [28] . It is a vector generalisa- 
tion of the Ito equation. 

Systems (1351) and (1361) seem to be new. They are vector generahsations of the 
corresponding scalar systems from the paper [21]. A coherent picture of such vector 
analogues of the Kac- Moody KdV-systems will be developed in a separate paper. It 
is interesting to note that the fifth order symmetry 



Ut = —9 Uxxxxx — SOm Uxxx — 45^^ U^x — (35iX^a; + 20u'^ + 5{U,U)) U^ 
-{Wuxxx + 20uux + 5{U,Ux))U 

of system (1551) is nothing but a vector generalization of the Kaup-Kuperschmidt 
equation. Indeed, if the vector part is absent (i.e. U = 0) then system ([35D becomes 
trivial and (1371) turns out to be the Kaup-Kuperschmidt equation. 

Our attempts to find vector analogues for the Sawada-Kotera equation among 
systems with one vector and one scalar unknowns were unsuccessful. It turns out 
that any non-triangular, non-degenerate system of fifth order with A = 2, having a 
symmetry of seventh order, is a symmetry of one of ([SI]), (151|) . ( 1551) or (I5B]) . 

Appendix 1. 

Definition of Jordan algebra: A vector space J equipped with a product o : 
J X J ^ J such that 



10{U, Uxxx) - iHUx, Uxx) - lOuxiU, U) 



20u{U,Ux), 



(37) 



XoY = YoX, o{Y oX) = {X^ oY)oX 



12 



for any X,Y G J is called Jordan algebra. 

If * is a multiplication in an associative algebra then XoY = X*Y + Y*Xis 
a Jordan product. 

Definition of Jordan triple system: A vector space J equipped with a triple 
product {■,■,■}: J X J X J ^ J such that 

{X,Y,Z} = {Z,Y,X} 

{X, Y, {V, W, Z}} - {V, W, {X, Y, Z}} = {{X, Y, V}, W, Z} - {V, {Y, X, W}, Z} 

for any V, W, X,Y, Z & J is called Jordan triple system. 

Any associative algebra with a multiplication * produces a Jordan triple system 
with respect to the triple product 

{X,Y,Z} = X *Y * Z + Z *Y * X. 

Appendix 2. 

For different ansatze of a vector equation and its symmetry, like equations fllSp 
and f|T6|) . the algebraic conditions for the unknown coefficients vary drastically in 
their size. For the nine combinations of 3 values A = 2, 1, ^ and 3 pairs of differ- 
ential order of the equation and its symmetry ((2+3), (3+5), (5+7)) the table below 
shows the number of unknowns, the number of conditions and the total number of 
terms in these conditions for the resulting nine over-determined systems of algebraic 
conditions. 



order 


2 + 3 


3 + 5 


5 + 7 









2 


8 


^ of unknowns 


A = 2 





4 


8 


# of equations 







4 


8 


total ^ of terms 




3 


10 


31 


^ of unknowns 


A= 1 


5 


26 


198 


# of equations 




9 


121 


3125 


total ^ of terms 




10 


33 


107 


^ of unknowns 


21 


129 


927 


# of equations 




80 


1603 


52677 


total 7^ of terms 



Table 1. A comparison of the size of the symmetry conditions for single vector 

equations. 
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